Lemaitre-Tolman-Bondi cosmological models, smoothness, and positivity of the 

central deceleration parameter 
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We argued in a previous paper [R. A. Vanderveld et al. 2006] that negative deceleration parameters 
at the center of symmetry in Lemaitre-Tolman-Bondi cosmological models can only occur if the 
model is not smooth at the origin. Here we demonstrate explicitly the connection between non- 
smoothness and the failure of positivity theorems for deceleration. We also address some confusion 
that has arisen in the literature and respond to some recent criticisms of our arguments. 



I. INTRODUCTION 

In a previous paper [1] (henceforth VFW) we studied 
spherically symmetric, dust cosmological models, which 
are described by the Lemaitre-Tolman-Bondi (LTB) met- 
ric. In particular we discussed the potential for such mod- 
els to reproduce the apparent acceleration of the Uni- 
verse, and how some models evade positivity theorems 
for the central deceleration parameter. Our results were 
criticized in the recent paper [2| by Krasinski, Hellaby, 
Celerier, and Bolejko (henceforth KHCB). Most of these 
criticisms are due to misinterpretations and are invalid. 
In this short note we clarify some of the issues involved 
and respond to the criticisms. 



II. TWO KEY MISCONCEPTIONS IN KHCB 
A. Misidentiflcation of deceleration parameter 

A key issue that arises is the definition of the de- 
celeration parameter. Throughout VFW we defined 
the deceleration parameter q(z) as a function of red- 
shift z to be what would be computed by observers 
from observed luminosity distances Dl(z) by assuming a 
spatially flat 1 Fricdmann- Robertson- Walker ("FRW" in 
VFW and "FLRW" in KHCB) cosmological model. This 
is discussed explicitly around our Eq. (2.10). Specializ- 
ing to z = gives the observed deceleration parameter 
at the location of the observer, go = <z(0). 

This discussion applies to spherically symmetric space- 
times, with the observer at the symmetry center. For 
more general inhomogeneous spacetimes, some prescrip- 
tion for angle averaging is required in order to define a 
deceleration parameter. One such definition at z = is 
given in Ref. [3j], which in spherical symmetry coincides 
with the definition of go given above. Slightly different 



1 If one instead allows for a spatially curved geometry, the deceler- 
ation q(z) at finite z can be altered, but the central deceleration 
q(0) at z = is unchanged, as we discussed in VFW. 



definitions, specialized to z = 0, are given in Hirata and 
Seljak Q (henceforth HS), namely the quantities which 
they denote by q% and q±. Those definitions also coincide 
with the above go in spherical symmetry. 

The starting point of our paper VFW was the following 
contradiction that existed then in the literature. Namely, 
the result of Ref. [3| specialized to spherical symmetry 
showed that go > always for any smooth LTB model. 
Similar results in HS related to their quantities q% and 
<?4, specifically the argument from their Eqs. (18) - (41), 
give the same result when specialized to spherical sym- 
metry. On the other hand, the explicit LTB solutions 
constructed in Ref. 0] and elsewhere (see VFW for ref- 
erences) have <?o < 0. We pointed out in VFW the res- 
olution of this apparent contradiction: the solutions of 
Ref. Q are not smooth at their center, and therefore vi- 
olate the smoothness assumptions used in deriving the 
positivity theorems of Refs. [3|, |4| . 

The first key error made in the comments of KHCB 
is a misidentiflcation of which definition of deceleration 
parameter is involved, wherein they focus on another def- 
inition given in HS - related to the derivative with respect 
to a fluid element's proper time of the expansion of the 
fluid - which HS denote by q% and KHCB denote by qHS- 
One of the main criticisms in KHCB is that we were con- 
fusing different definitions of q in VFW, using both go 
defined above and q\, but this is not the case as qi was 
never referred to in VFW. In addition, when we cited 
HS we explicitly stated which definition of acceleration 
parameter was involved, saying "the local expansions of 
. . . Hirata and Seljak show that go is constrained to be 
positive." We then proceeded to specify the definition of 
deceleration parameter go [our Eq. (2.16)] and reviewed 
the corresponding positivity theorem. 

Having focused on the incorrect deceleration param- 
eter qi, KHCB then argued, correctly, that there is no 
contradiction between the positivity theorem of HS for 
qi and the explicit models of Ref. [5[. This is because, 
although the definitions of go and gi coincide for smooth, 
spherically symmetric models at z — 0, they do not coin- 
cide for models which are not smooth, and in particular 
the models of Ref. @ have gi > even though g < 0. 
While KHCB are correct about this point, it is not rele- 
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vant to a discussion of VFW. 

We emphasize that we used only a single definition 
of deceleration parameter throughout our paper, and 
KHCB's claim that we "intermittently used the same 
name. . .for two different quantities" is incorrect. Al- 
though KHCB give two different formulae that they claim 
we use, their Eqs. (2.12) and (2.13), in fact the first for- 
mula is just the second formula evaluated at z = 0. 

There are further criticisms in Sec. VI of KHCB re- 
lated to our computation of the deceleration parameter. 
After their Eq. (6.1) KHCB say that "this makes it clear 
that the shear and vorticity terms . . . are also evaluated 
at z = 0, though this is not clear in either paper." In 
fact this is explicitly stated in our Ref. [H, where the 
first sentence in the paragraph containing Eq. (3) reads 
". . .in terms of the density, 4- velocity and velocity gradi- 
ents of the cosmological fluid evaluated at the observer's 
location V." It is also stated explicitly in VFW, where 
qo is defined to be the "central deceleration parameter" 
before Eq. (2.16), so it is clear that everything is evalu- 
ated where the observer lives, at z = 0. The fact that the 
gradients of the velocity are evaluated at z — also in- 
validates KHCB's criticism in the first paragraph of their 
Sec. VID. 



B. Skepticism about relevance of differentiability 
of solution 

Because there is no contradiction involving the deceler- 
ation parameter q\ , KHCB believed that our explanation 
in terms of non-smoothness was invalid. They did appre- 
ciate that a contradiction still existed for the decelera- 
tion parameter q$, which they summarized in their Sec. 
VIA. However, they discounted our resolution of the con- 
tradiction, and instead argued that the formula for the 
deceleration parameter derived in the earlier paper Q 
was "erroneous or erroneously interpreted in the VFW 
paper." 

In order to address this skepticism, we give here a de- 
tailed and explicit demonstration of the connection be- 
tween smoothness of the geometry and the positivity the- 
orems of Refs. [1, H[. Although such a demonstration is 
not strictly necessary since the argument of VFW is com- 
plete and self contained, it may nonetheless help to dispel 
confusion about this issue. 

A basic ingredient of those positivity theorems is a 
local covariant Taylor expansion of the fluid's 4-velocity 
about a given point x a , of the form 



Ua'(x') — g aa ' (x, x') \u a {x) + u a/3 (x)a- a (x, x')a-p{x, x') 

(2.1) 



+ -u at3j (x)a- a (x, x')<7 ] p{x, x')a-, 1 (x, x') 



where 



.a/3 



(x) = -V<- a ufl(x) and 



a/3~f 



(X) 



V^ a V /3 u 7 ^ ) (x). Here g aa '(x,x') is the parallel transport 
bivector and a(x,x') is Synge's worldfunction; see Ref. 



Q for details. Clearly, in order for this expansion to be 
valid, the coefficients u a , u a ", and in the expansion 
must exist, and in particular the symmetrized second co- 
variant derivative of the 4-vclocity 



V( ct V /3 u 7) 



(2.2) 



evaluated at x a must exist. We now show that this quan- 
tity fails to exist at the center of symmetry for the LTB 
models with qo < 0, which explains why the positivity 
theorems are inapplicable to such models. 

The metric of LTB models can be written in the general 
form 

ds 2 = -dt 2 + e 2a ^dr 2 + e 2/3 ^r 2 (d9 2 + sin 2 9dip 2 ), 

for some functions a(t, r) and 0(t, r). The fluid 4-velocity 
is u a — — {d£) a . We assume the differential structure on 
the manifold associated with the coordinates (t,x,y, z), 
where (x, y, z) are given in terms of (r, 6, ip) by the usual 
formulae for spherical polar coordinates 2 . Then, it can 
be seen that both u a and u a are smooth tensor fields. 
However the metric g a p and connection V Q need not be 
smooth, so covariant derivatives of u a need not exist. 

A convenient piece of the second derivative to look at 
is h a = (g a P +u a u l3 )g' tS V(0V 1 us). From the metric > 
we then obtain 



h a dx a = 



dr, (2.4) 



where dots denote derivatives with respect to t and 
primes denote derivatives with respect to r. We now 
expand the functions a and f3 as 

o(i,r) = a (t) + a 1 (t)r + a 2 (t)r 2 + 0(r 3 ), (2.5) 
(3(t,r) = /3 (t)+/3i(t)r + /3 2 (<)r 2 + O(r 3 ). (2.6) 

For all LTB models we have ao = /?o [see Eq. (2.20) of 
VFW] , and so we get near r — that 



h a dx c 



7di(t)-2/3i(i)+0(r) 



dr. 



(2.7) 



If the quantity in square brackets is nonzero, then h a 
is nonzero in the limit r — > 0, and points in the radial 
direction. Therefore it has a direction dependent limit as 
r — ► 0, i.e., the limit does not exist. (Correspondingly, 
higher order derivatives of this quantity diverge.) 



2 For this differential structure the four velocity is smooth (C°°) 
while the metric is C 1 but not C 2 . It is possible to find an alter- 
native differential structure [using a coordinate transformation 
of the form f = r + f(t)r' 2 + . . . for some function f(t)] for which 
the metric is C 2 but not C 3 and the four velocity is C 1 but not 
C 2 . However there is no choice of differential structure for which 
the metric is smooth, as can be seen by computing coordinate in- 
variants like VaV a i?, where R is the Ricci scalar, which diverge 
as r -* 
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The coefficient in Eq. (|2.T[) will generically be nonva- 
nishing for the type of models in Ref. [l[ with q < 0. 
Such models were characterized by having non- vanishing 
linear terms in the expansions of the bang time function 
to (r) = too + toiT + t^r 2 + 0(r 3 ) and curvature function 
k(r) = ko + k\r + k2r 2 +0(r 3 ). Smoothness requires that 
k% = ioi = 0. For example, for the case k(r) = we have 



h n dx a 



3^ 



0(r) 



dr. 



(2.8) 



Therefore we see explicitly that for the models with go < 
the second covariant derivative of the 4-velocity does 
not exist, explaining why the positivity theorem does not 
apply. 



III. CRITICISMS RELATED TO CHOICES OF 
TERMINOLOGY 

Many of the criticisms in KHCB are not directed at 
our mathematical results, but instead are related to the 
choice of terminology we employ to describe our results. 
In most cases the criticisms are based on misinterpre- 
tations and are unfounded, but some of their criticism 
has validity, since some of the terminology we used was 
poorly chosen and apt to lead to confusion. We now dis- 
cuss and clarify the relevant choices of terminology: 

1. "Singularity:" In VFW we discussed in detail the 
nature and implications of the lack of smoothness 
at the origin of the LTB models with qo < 0. 
We called these locations singularities, in the loose 
sense that there were some coordinate-invariant 
quantities which become singular there. This was 
a poor choice of terminology, since "singularity" 
is usually used in the general relativity literature 
to mean geodesic incompleteness, which does not 
apply here, and/or a divergence of the Riemann 
tensor, which also does not apply here. Although 
we did not claim that the singularity was a curva- 
ture singularity, our terminology was confusing, as 
correctly pointed out by KHCB. 

KHCB also object to our describing the singularity 
as "weak." While it is true that there are variety of 
different terminologies for classifying singularities 
in use in general relativity, "weak" and "strong" 
are now in common use as referring to whether par- 
allel propagated components of the tidal distortion 
tensor diverge or not. The precise definitions of 
weak and strong are given, for example, in Ref. Q . 
The spatial origin in the non-smooth LTB models is 
weak in this sense (rather trivially, since the compo- 
nents of the Riemann tensor are finite). Therefore 
our description is appropriate and the criticisms of 
KHCB on this point are unfounded. 

2. "Unphysical:" In VFW we asserted that the non- 
smooth LTB models are unphysical. KHCB dis- 



agreed by arguing that one can smooth out the cen- 
tral singularity easily at z < z s for some small red- 
shift z s without changing the predictions at larger 
z. We agree with this point, and in fact we men- 
tioned it in the concluding section of VFW. The 
possibility of performing such smoothing is well- 
known, see, for example, the numerical studies in 
Ref. [7|. However, the non-smooth models are still 
unphysical when the observer is placed at the cen- 
ter, in the following sense: Such models have con- 
siderably more fine tuning than smooth models, 
since it is unlikely for an observer to live exactly 
on such a point in the density distribution, and 
especially since to smooth out the singularity in a 
manner compatible with luminosity distance obser- 
vations requires the introduction of a new, artificial, 
very small lengthscale. 

3. "Inverse Problem:" In Sec. II of VFW we discussed 
the straightforward computational procedure for 
obtaining the luminosity distance D^(z) from the 
bang time function to(r) of a zero-energy, LTB 
model. In Sec. Ill we discussed the what we called 
the "inverse problem," by which we meant simply 
the inverse process of attempting to find to(r) from 
a specified Dl(z). KHCB appear to interpret the 
phrase "inverse problem" in a different and much 
more general sense, and therefore their criticism on 
this point is unfounded. 

4. "Effective Equation of State:" In VFW we defined 
the effective equation of state parameter w e g(z) to 
be the equation of state that is obtained from the 
data when assuming the usual framework of a flat 
FRW cosmology. We by no means said that it is 
the equation of state underlying the LTB model 
itself (which is of course that of pressureless dust). 
KHCB seem to be confused about this point in their 
Sec. VIB, and appear to believe that we intended 
a literal interpretation of w e g(z). 



IV. FINDING LTB MODELS THAT YIELD A 
SPECIFIED LUMINOSITY DISTANCE AS A 
FUNCTION OF REDSHIFT 

In Sec. Ill of VFW, we explored the problem of trying 
to find an LTB model that would have the same lumi- 
nosity distance Dl(z) as a given FRW model with or 
without dark energy. We called this the "inverse prob- 
lem." We found this problem to be complicated by a 
generic critical point in the differential equations to be 
solved. However we were still able to show that "trans- 
critical" LTB models could be constructed. In Sec. V of 
their paper, KHCB criticize our discussion on a number 
of issues. Their criticism is largely unfounded. 

First, KHCB appear to have misread our paper on one 
important point. In Sec. V, they say "VFW . . . argue 
that only [dust] FRW models have ... [a] critical point," 
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and "[VFW] say that [dust] FRW models provide ex- 
amples of transcritical solutions . . . but they fail to find 
other viable [transcritical] solutions." These two asser- 
tions are incorrect, since in our Sec. IIIB we explicitly 
constructed a variety of transcritical solutions that corre- 
spond to choices of Dl (z) other than dust FRW. The ori- 
gin of KHCB's misinterpretation might be the fact that 
we gave the explicit formula angular diameter distance 
for dust FRW models after Eq. (3.9). However, this was 
given only as an illustrative example. 

Second, KHCB appear to misinterpret our statements 
about the generality of critical points. They argue that 
all light rays in physically reasonable cosmological mod- 
els do encounter critical points (local maxima of angular 
diameter distance as a function of redshift). We agree 
with this assertion, and we never claimed otherwise in 
VFW. After our Eq. (3.13) we say "only the special class 
of transcritical solutions can extend to infinite redshift." 
Although we did say these solutions were "special," we 
intended only to mean special within the mathematical 
class of solutions to our differential equations (3.4) and 
(3.5), and not within the space of physically acceptable 
LTB models. KHCB appear to have misinterpreted us 
on this point. 

Third, in VFW we argued that not every angular di- 
ameter distance function Da{z) can be realized by zero- 
energy-function LTB models. KHCB argued that generic 
choices of Da (z) are realizable. However, the evidence 
in favor of this assertion they presented is in the context 
of fully general LTB models, in which one varies both 
the energy function and the bang time function, a con- 
text different from that of VFW. For the context studied 
in VFW, where one varies only the bang time function, 
it is straightforward to explicitly confirm that not every 
Da (z) can be realized by studying LTB models that are 
linear perturbations of dust FRW and using the test out- 
lined after Eq. (3.25) of VFW. In the notation of VFW 
we choose the angular diameter distance to be 



where e < 1 is a small parameter and S(z) is an ar- 
bitrary function. We choose 8{z) to vanish in a neigh- 
borhood of z — and in a neighborhood of the critical 
point at z = 5/4. Then the linearized version of the 
differential equation (3.15) is V{(z) — j(z)S'(z), where 
we have parameterized the function V(z) as V(z) — 

3 - 2a/ITz + e{\ + z) 2 / 3 [3-2\/rTi] _1/3 Vi(z) and 
the function 7(2) is 



j(z) = — 



2^/3-2V7TT 7l (z) 



r-FRw(^) 



[l + eS(z)}, 



(4.1) 



(z + l) 2 /3 (4z (V7TT - 4) + 25 (Vz~+T - 1)) 

(4.2) 

with 71(2) = -2z 2 +(9V^TT- 19) z + 20 (y/z + T-l). 
The boundary conditions for the differential equation for 
Vi(z) are Vi(Q) = Vi(5/4) = 0, and it follows that a 
transcritical solution will be possible for a given choice 
of 8{z) only if J^ 4 r y(z)8'{z) = 0. Therefore there are 
many choices of angular diameter distance that are not 
realizable. 



V. CONCLUSIONS 

Because of recent criticism 0, and to prevent future 
misunderstandings, we have clarified some of the points 
made in VFW PJ . In doing so we have shown that many 
of the criticisms are unfounded. 
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